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We examine the effects of pure additive noise on spatially extended systems with quadratic 
nonlinearities. We develop a general multiscale theory for such systems and apply it to 
the Kuramoto-Sivashinsky equation as a case study. We first focus on a regime close to 
the instability onset (primary bifurcation), where the system can be described by a single 
dominant mode. We show analytically that the resulting noise in the equation describing 
the amplitude of the dominant mode largely depends on the nature of the stochastic 
forcing. For a highly degenerate noise, in the sense that it is acting on the first stable 
mode only, the amplitude equation is dominated by a pure multiplicative noise, which in 
turn induces the dominant mode to undergo several critical state transitions and complex 
phenomena, including intermittency and stabilisation, as the noise strength is increased. 
The intermittent behaviour is characterised by a power-law probability density and the 
corresponding critical exponent is calculated rigorously by making use of the first-passage 
properties of the amplitude equation. On the other hand, when the noise is acting on the 
whole subspace of stable modes, the multiplicative noise is corrected by an additive-like 
term, with the eventual loss of any stabilised state. We also show that the stochastic 
forcing has no effect on the dominant mode dynamics when it is acting on the second 
stable mode. Finally, in a regime which is relatively far from the instability onset, so that 
there are two unstable modes, we observe numerically that when the noise is acting on the 
first stable mode, both dominant modes show noise-induced complex phenomena similar 
to the single-mode case. 

Key Words: Stochastic partial differential equations, critical phenomena, multiple scale meth- 
ods, singular perturbations. 



1 Introduction 

External or internal random fluctuations are ubiquitous in many physical systems and can 
play a key role in their dynamics often inducing a wide variety of complex spatio-temporal 
phenomena. Examples can be found in several fields: From biology (such as stochastic 
resonance accounting for the sensitivity that certain species have, to weak but coherent 
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signals in noisy environments) and technological applications (e.g. high-temperature su- 
perconductivity [2 [51]), to fluid dynamics (e.g. Rayleigh-Benard convection commonly 
used as a prototype to study instabilities in systems out of equilibrium [7] , and waves in 
free-surface thin film flows [23]). Many of these natural phenomena and technological ap- 
plications can be described by model noisy spatially extended systems (SES), i.e. infinite- 
dimensional dynamical systems described through stochastic partial differential equations 
(SPDEs) with space-time dependence [H], such as the noisy Swift-Hohenberg equation, or 
the noisy Kuramoto-Sivashinsky (KS) equation. The spatio-temporal dynamics of noisy 
SES can be dominated by many curious phenomena, such as noise-induced spatial pat- 
terns [13] and noise-induced phase transitions [16] , Clearly, characterising the influence 
of noise in SES is crucial for the understanding of the inception and long-time complex 
spatio-temporal dynamics of physical systems, as well as for the control and optimisa- 
tion of technological processes. The identification and understanding of different regimes 
in parameter space, including the emergence of underlying scaling laws is of particular 
interest in analytical and computational studies of such systems. 

One of the most widely studied complex phenomena induced by noise effects is the 
transition between different observed system states as the noise strength is continuously 
increased beyond a critical value. A related question concerns the mechanisms by which 
fluctuations in the system interact with nonlinearities to induce ordered or regularised 
states (see [21 [H] ) . It is now generally accepted that one of the main mechanisms required 
to induce phase transitions (e.g. between ordered-disordered phases), is the presence of 
pure multiplicative nois^], even though the presence of combined multiplicative and 
additive noise has also been shown to induce phase transitions 26 5 [54]. On the other 
hand, complex phenomena such as stabilisation effects (but not phase changes) can be 
induced by pure additive noise as reported in recent numerical investigations [1 71 1181 134] . 

In this study, we investigate the effects of pure additive noise on unstable SES that are 
close to the primary bifurcation ("instability onset"). A first stab at this problem was our 
recent study in [41] which focused on the influence of highly degenerate noise, a particular 
type of additive noise, on the KS equation - the noisy KS equation is a paradigmatic 
model for a wide spectrum of physical settings. In [41^ it was shown that being close to 
the instability onset allows for a global description of the KS equation in terms of a single 
dominant mode. In the present work, by means of a multiple scale analysis for general SES 
with quadratic nonlinearities and by appropriately extending the singular perturbation 
methodology in |38j for such SES, we obtain an amplitude equation for the dominant 
mode which then enables us to describe analytically the behaviour of this mode and 
explore systematically the effects of noise. As an application of the general methodology 
we develop we examine the noisy KS equation. We offer a derivation of this equation for 
a hydrodynamic system, that of a thin-liquid film flowing over a topographical substrate 
obtained asymptotically from governing equations (Navier- Stokes) and wall and free- 
surface boundary conditions. In other contexts, the noisy KS equation has been adopted, 

Multiplicative noise enters the system through external effects, i.e. via the presence of 
noise on a controlling parameter or through noisy fluctuations in the boundary conditions; on 
the other hand, additive noise arises from fluctuations of an internal origin such as thermal 
fluctuations, which although they are present at the nanoscale, they can often have macroscopic 
effects (e.g. [33]). 
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for instance, as a model for sputtering processes used to produce thin films in materials 
applications including nanostructuring solid surfaces using ion beam erosion - in addition, 
a possible use of noise as a control tool is also suggested and explored computationally, 
see P [HI EH [291 [30] . 

We observe that the dynamics, and in particular the statistical properties of the fluc- 
tuating dominant mode, are largely dependent on the nature of the applied noise and 
its strength. More precisely, we consider a degenerate noise, in the sense that it is acting 
on the subspace of stable modes alone, and scrutinise its effects on the dynamics of the 
unstable mode. Several situations arise depending on what stable mode (or modes) the 
noise is acting on. For example, when the noise is acting on the first stable mode only, 
the governing equation for the amplitude of the dominant (unstable) mode is reduced 
to a purely temporal Stuart-Landau (SL) model with pure multiplicative noise. The nu- 
merical simulations presented in our previous study [41] suggested that this type of noise 
is in turn responsible for the different critical transitions that the dominant mode un- 
dergoes as the noise strength is increased, including an initial state of finite fluctuations, 
an intermittent on-off state characterised by a burst-like dynamics, and a completely 
stabilised state. It is noteworthy that such critical transitions exhibit universal underly- 
ing scaling laws, in the sense that the observed critical exponents are also ubiquitously 
found in many different physical systems. It is also noteworthy that such observed on-off 
intermittency is actually a crucial feature of many nonlinear systems close to criticality, 
and reflects a transition from order/coherence to a disordered state. 

Here, we rigorously derive the statistical properties of such intermittency in terms 
of a universal law that does not depend on the particular model under consideration. 
We further offer new results and insights on the influence of noise on SES compared to 
our previous study in [41) . More specifically, when the noise acts on the whole set of 
stable modes, we show that the multiplicative term in the SL equation is corrected by 
an additive-like term, and any stabilisation effect is eventually lost. Additionally, we find 
that the noise can be filtered out by the nonlinearities when it is acting on the second 
stable mode alone. In this case, the dynamics of the dominant mode is not affected by 
the stochastic forcing applied to the system. Our analytical findings are in full agreement 
with numerical experiments of the KS equation, thus providing a complete picture of 
the relevant effects of additive noise on SES with quadratic nonlinearities close to the 
instability onset. 

Following the development of analytical and computational descriptions of the phe- 
nomena near the primary bifurcation point, it is important to consider the dynamics 
beyond criticality where the multiscale analysis is no longer valid. To achieve this we 
revert to numerical experiments and compute the resulting spatio-temporal dynamics 
when two linearly unstable modes are active. Surprisingly, when we consider a highly 
degenerate noise acting on the first stable mode alone, the dynamics of the two domi- 
nant modes are found to follow the same critical transitions as with the case close to the 
primary bifurcation where only one unstable mode is present. 

In f}2]we present the theoretical framework for noisy SES close to the primary bifurca- 
tion. In EI2.1l we derive the generic amplitude equations for the unstable mode components 
using multiscale techniques. Our theory is then applied to the noisy KS equation in <J3] 
and both analytical and numerical results for the different cases of noise considered are 
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presented in ^3 . 1 1 where we take the noise acting on the first stable mode only, in ^3.21 
where the noise acts on the first and second stable modes, and in £)3.3[ where the noise 
acts on the second stable mode only. A numerical study in a regime relatively far from 
the instability onset is presented in Sj4] Finally, our results are summarised in <j5l 

2 Theoretical framework: Noisy SES close to the primary bifurcation 

Typically, noisy SES can be described through SPDEs of the following generic form: 

d t u = Cu + 7"(it, Vji, V 2 u, . . .) + cr£(r, t), (2.1) 

where C is usually a linear differential operator with constant coefficients, and J 7 is a 
nonlinear function of its arguments. The field u(r, t) describes the magnitude of a quantity 
of interest in the system, and we also include the presence of a stochastic additive forcing, 
given by the random variable £ [r, t), with a its strength. The complexity of SES of the 
form (|2.ip and their dynamic^ is such that it is often quite difficult, if not impossible, 
to analyse them directly either analytically or even numerically. It is therefore desirable 
to obtain an appropriate low-dimensional description representing the dynamics at a 
coarse-level and which can capture most, if not all of the essential dynamic features of 
the original high-dimensional model describing the actual problem. 

Quite often in physical systems, the dominant nonlincarity is a quadratic one whose 
functional form is dictated by simple symmetry considerations. For example, in the con- 
text of free-surface thin-film flows, the dominant nonlinearity is ud x u, associated with 
the interfacial kinematics due to mean flow [45] (the only other dominant nonlinearity, 
u 2 , is easily ruled out for systems whose spatial average is conserved, i.e. dt(u) x — 0). We 
then consider SES with a quadratic nonlinearity, random forcing and a spatially uniform 
rest state. The spatial domain is taken to be [-L, L] d where d + 1 is the number of space 
dimensions. We note that even though the multiscale techniques we shall be developing 
in this study could, in principle be extended to multidimensional systems, this extension 
is non-trivial and beyond the scope of the present study (in fact, the theory of SPDEs 
in dimensions higher than one is currently not well developed, at least for space-time 
white noise). For this reason, we will consider Eq. (12.11) in one spatial dimension (d = 1). 
We also assume that zero is an eigenvalue of the linear operator C with the correspond- 
ing eigenspace being finite dimensional, and the remaining eigenvalues being in the left 
complex half-plane. 

Let us introduce a bifurcation parameter e, such that for e = the system approaches 
its rest state as t — > oo. For e > 0, a bifurcation occurs leading to a finite number of modes 
that are linearly unstable. Therefore, we assume that F{u, Vu, V 2 u, ■■■) = e 2 Ju+B{u, u), 
where J and B(u, u) are a linear differential operator with constant coefficients and a 
quadratic nonlinearity (a symmetric bilinear map), respectively, so that for e«l, there 
appears at least one eigenvalue of the perturbed operator C + e 2 ^ in the right half-plane. 

2 It should be emphasised, that already in the absence of noise (a = 0), the presence of a wide 
range of spatial and temporal scales coupled with each other through nonlinearities may give 
rise to a complex dynamics such as transitions between different patterns and spatio-temporal 
chaos [7]- 
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Equation (|2.1[) is then written in the form: 

d t u = Cu + e 2 Ju + B(u, u) + eer£(x, t). (2.2) 

For the particular distinguished limit that we will consider in this paper, and which 
leads to the most interesting coarse-grained dynamics, we have rescaled the noise term 
cr£ so that its strength is of O(e), in comparison to the 0(e 2 ) distance from criticality. 
Different distinguished limits, relevant when different assumptions on the nonlinearity 
and the structure of the noise are made, can also be considered, e.g. [3l [6]. The field 
u can then be projected onto the set of eigenfunctions {ek(x)} for k = 1, . . . , oo of the 
linear operator C, i.e. u(x,t) — ^2kUk(t)ek{x). For < e <C 1, the system is close to the 
bifurcation point and is described by the presence of a finite number of unstable modes, 
the "dominant modes" . In this regime, Eq. (12.21) has two widely separated time scales, 
corresponding to the (stable) fast and (unstable) slow modes, allowing us to derive an 
amplitude equation for the slow-dynamics of the dominant modes only, which belong to 
the null space Af of the linear operator C. 



2.1 Amplitude equation reduction 

We start with the generic stochastic equation (|2.2[) on the domain [-L, L]. For e -C 1, the 
number of linearly unstable modes is given by the dimensionality of Af. We are interested 
in the dynamics of the dominant modes when the stable modes, «fc(t)efc for e^ E Af , 
are randomly forced. Here, Af ± represents the subspace of fast modes orthogonal to Af. 
The noise term in Eq. (|2.2I) is therefore written as: 

£(x,t) = Qk/3 k {t)e k (x), (2.3) 

where the index k represents the modes such that e& <G Af ■ The variable $k(t) in the 
above equation corresponds to uncorrelated white noise, (/3 m (t)/3 n (t')) = 5 mn 5(t—t'), and 
qk represents the wave number dependence of the noise. Considering now the behaviour 
of small solutions at time scales of 0(e -2 ), we define u(x,t) — ev(x,e 2 t) and use the 
scaling properties of the white noise to transform Eq. (|2.2[) to: 

d t v = e- 2 Cv + Jv + e^B(v, v) + e^a^x, t). (2.4) 

A detailed rigorous derivation of the pertinent amplitude equation for the dominant mode 
in the case where J is the identity operator can be found in Ref. [4] . We shall extend this 
previous formalism here for the general case of J being any linear (differential) operator 
that commutes with £, i.e. both operators have the same eigenfunctions, and for a finite 
dimensional kernel, i.e. dim(A/") = A^o- 

To obtain the amplitude equation for the dominant modes we consider a finite dimen- 
sional truncation of the above system and keep M fast modes in the series expansion, 
so that the total number of modes in the expansion is N = No + M. We also consider 
Eq. (|2.4j) in a bounded domain, e.g. [— tt, tt], and we project the field v onto AT to get 
Vi = V c v, where V c is the corresponding projector to the null space, and onto its or- 
thogonal subspace Af 1 ' to get v± = V s v, where V s = 1- — V c with I being the identity 
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operator. The resulting system of equations reads: 

dtvi =^i + 2£- 1 P c B(« 1 ,« i ) + e- 1 ? c 8( ! ; i , t .i), (2.5 a) 

d t v x = (J- e- 2 C)v ± + e- x VJ5{v, v) + e _1 o£(»> *)• (2-5 6) 

By applying now techniques from homogenisation theory |38j and assuming the condition 
for the bilinear map of VJB( e k> e fe) = for k = 1, . . . , N, we can derive the homogenised 
(i.e. amplitude) equations that describe the dynamics of the unstable modes. The detailed 
derivation for a multidimensional null space is given in Appendix |A"1 

In the case of ID null space, we dehne the amplitude A for the dominant mode as Vi = 
A(t)ei, and its corresponding equation is given by the SL equation with multiplicative 
noise (see Appendix lAl for details): 

A = (ji + no 2 ) A - l2 A z + a^ la a 2 + lm A 2 W(t), (2.6) 

with W{t) is a white noise, and where the different coefficients are given in Eq. (|A 21|) 
from Appendix [S] We therefore see that the coefficients explicitly depend on the nature 
of the noise, i.e., the values of q^, and the properties of the quadratic nonlinearity, which 
are dictated by B n ^ given in Eq. (|A3|) . 

In the case of a two-dimensional (2D) null space, the solution is expanded as v\ — 
a\{t)ei + a2(t)e2, and the coupled amplitude equations we obtain are of the form: 

di = (ji + 2 7 icr 2 )ai - j 2 aiA 2 + 2a^j a a 2 + j m A 2 Wi(t), (2.7 a) 

d 2 = (, n + 2 7l( x 2 )a 2 - l2 a 2 A 2 + 2o^J la o 2 + lm A 2 W 2 (t), (2.7 b) 

where we have defined the global amplitude A(t) = \/a\ + a\, and the coefficients 71, 
72, 7a, and 7 m are given by Eqs. (IA21I) . 



2.2 Analysis of the amplitude equation 

The stationary statistical properties of the amplitude dynamics can be studied by solving 
the corresponding stationary Fokker-Planck equation. In the ID case, the stationary 
probability density function (PDF) for the random variable A for natural boundary 
conditions is given by [THl [21] : 



N c f A v(z) 



P (^) = T7^exp/ -^dz, (2.8) 



where v(A) — (j\ + r yio 2 )A — 72 A 3 and g(A) = ay/ j a o 2 + 1mA 2 , yielding in our case 
P(A) = N C {A 2 + X 2 a 2 ) ai/2 exp (~^A 2 ), (2.9) 

with: 

, > 1+ (71 +72A 2 )ct 2 72 , 01 „s 

a 1 [a) = 2 1, n(p) = o> ( 2 - 10 ) 

•y m (r -fni<J 

where we have defined the parameter A = 7 Q /7m, and the normalisation constant N c 

is given as A^ c = 2//^ /r A (ai), where T A (ai) = J^x^^ix + ^X 2 a 2 ) a ^ 2 e~ x Ax. It 

should be noted that the statistical properties of the fluctuating dominant mode depend 
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on both the strength of the noise a and the values of qk (see $3] below for the example of 
the KS equation). 

For a 2D null space, the stationary joint PDF for the two variables, G(ai,a,2), can 
also similarly be obtained by computing the corresponding stationary 2D Fokker-Planck 
equation. This yields to: 

G(a 1 ,a 2 ) oc (a 2 + a 2 + A 2 ) Q ' l/2 exp [-//(a? + a%)], 

where the parameters a[ and fj,' are obtained from the expressions in Eq. (|2.10[) by 
replacing 71, j a , and j m with 271, 47 a , and i-j m , respectively. The interesting point 
now is to study the behaviour of the PDF, P(A), corresponding to the global amplitude 
A = \/ a\ + a\. To this end, we move to a polar coordinate system (A, 9) by applying the 
transformation y\ = Asm9 and Z\ — Acos9, and we impose the normalisation condition 
between both distributions, G{a\, a2)daida2 = P(A, 9)dAd9, giving rise to a function of 
the form: 

P(A) cx A(A 2 + \ 2 )<' 2 exp (-^1'A 2 ). (2.11) 
Let us now apply the theory outlined above on a prototype model equation. 

3 Case study: the noisy Kuramoto-Sivashinsky equation 

Consider the noisy KS equation: 

d t u — — (d 2 + vd x )u — ud x u + er£, (3.1) 

normalised to 27r-periodic domains so that < v = (tt/L) 2 , where 2L is the original 
length of the system. The equation corresponds to an important class of SES, active- 
dissipative nonlinear media, whose main features are the presence of mechanisms for 
instability/energy production at long scales (d 2 u) and stability/energy dissipation at 
short scales {d x u). Both without and with the noise term, the KS equation has attracted 
a lot of attention since it appears in a wide variety of applications and physical phenom- 
ena. These include instabilities of flame fronts, models of collisional trapped ion modes 
in plasmas, dissipative turbulence, interfacial instabilities in two-phase flows, reaction- 
diffusion systems, the control of surface roughness in the growth of thin solid films by 
sputtering, step dynamics in epitaxy, the growth of amorphous films, and models in pop- 
ulation dynamics [HI HH ESI EHl SZ1 SS] ■ A detailed derivation of the noisy KS equation 
in the context of thin- film hydrodynamics is given in Appendix [Bj 

The KS equation represents an ideal candidate for our studies due to the extensive 
rigorous and computational results available. It is often used as a paradigmatic par- 
tial differential equation (PDE) that has low dimensional behaviour producing complex 
dynamics such as chaos - [HI [201 1221 133 EH1 [S3] • The existence and uniqueness of solu- 
tions for the stochastic KS equation have been proven in [TUj . The effect of weak additive 
white noise on transitions between stable fixed-point solutions and stable travelling-wave 
solutions of the KS equation has been considered numerically in [50) and the relationship 
between noise-induced transitions and the underlying attractors are explored. 

By assuming v = 1 — e 2 and a = e<r, Eq. (|3.ip is rewritten as: 

d t u = -{dl + dt)u + e 2 dtu - ud x u + e<j£. (3.2) 
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We can therefore read off all the different terms in Eq. (|2 .2[) as C = — <9 2 — d^, J — 
d*, and B(u,u) = ~ud x u. It is important to emphasise that in this setting, the term 
e 2 d^u represents a linear instability term that can destabilize the dominant modes of 
the equation. Note that it is controlled by the parameter e 2 which measures the distance 
from bifurcation. In this sense, a decrease in v = 1 — e 2 below the bifurcation point v = 1 
(and hence increasing e > 0) reduces the linearly stabilizing term d^u and therefore 
destabilizes the dominant modefjl. In our analysis we assume solutions of zero mean and 
we shall consider both the case of Dirichlct boundary conditions (DBC), i.e. u(—-K,t) = 
u{ir, t) = 0, and periodic boundary conditions (PBC). In the former case, the null space 
of C is ID: 

N{C) = span{sin(-)} , 

and the solution can be expanded in the basis {efc(x) = Ck sin (kqx)}, where q — tt/ L = 1 
and Cfe's are normalisation constants. On the other hand, for PBC the null space of C is 
2D: 

Af(£) = span{sin(-), cos(-)} , 

and the solution is expanded in the exponential Fourier basis {ei-(x) = Ck exp (ikqx)}, 
for fc = 0,±1,±2, . . .. By applying the multiscale formalism developed in we then 
obtain that the amplitude equation for the dominant mode is given either by Eq. (|2.6| 
for DBC or by Eqs. (|2.7p for PBC, with j\ = 1 and the following coefficients: 



-,2 
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Q2 , 1 9 » ( n ~ 1 n+l\ 



8 A2(A2 + A3) 8 A„ \A„ + A n _i A n + A Jl+ 



72 = 4' 7m = ^' {3 - 3&) 

_ 1 gfgf ly> 2^ , Qn-l \ /„ o \ 

' a " 8 A 3 (A 2 + A3) 2 + 8 ^ n \ A„ +1 (A„ + A n+1 ) 2 + A„_ 1 (A„ + A„_i) 2 ^ ° ] 

where Xk = k A — k 2 . From the above expressions it is interesting to note that several 
situations arise for different cases of noise. For example, for a highly degenerate noise 
acting on the first stable mode only (Case I, qu = ^.2), there is only multiplicative noise 
(7 a = 0). On the other hand, when the noise acts on both the second and third mode 
(Case II, qk — &k.2 + ^,3), the multiplicative noise is corrected by an additive-like term 
(7a 7^ 0). Finally, when the noise acts on the third mode only (Case III, q^ = Sk.3), the 
resulting amplitude equation contains no noise ("/ rn — 7 a =0). 

In the following, we analyse analytically and numerically the dynamics of the amplitude 
of the dominant mode A{t) = (t) | by considering the above three different cases of noise 
separately. To solve numerically the KS equation we adopt a pseudo-spectral method 
for the spatial derivatives that uses the Fast Fourier Transform (FFT) to transform the 
solution to Fourier space. The nonlinear terms are evaluated in real space and transformed 
back to Fourier space by using the inverse FFT. The solution is then propagated in time 

3 Note that when the KS equation is written in this form, multiscale analysis can be used in 
order to obtain the amplitude equation in a rigorous and systematic way. See appendix A for 
details. 
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Figure 1. Typical spatio-temporal evolution of the noisy KS equation (|3.1j) with PBC for 
a = 10 and e = 0.025, and imposing a highly degenerate noise acting on the first stable 
mode (Case I). 

by making use of a fourth-order Runge-Kutta scheme. In our simulations we have chosen 
a time step of At = 0.1, and the numerical error has been controlled by monitoring some 
constant quantity, namely the spatially averaged solution (u) = ^ /_£ u{x, t)dx, so that 
with the chosen time step it remains constant during all the numerical experiments. 

3.1 Case I. Amplitude equation with multiplicative noise: Noise induced 

critical state transitions 

We start by considering the case of highly degenerate noise acting on the first stable 
mode only so that the resulting amplitude equation for the dominant mode contains 
a pure multiplicative noise term (~/ a — 0). Typical snapshots of the spatio-temporal 
evolution of the KS Eq. fl3.1[) in this noise setting and subject to PBC with a — 10 and 
e = 0.025 are depicted in Fig. Q] The interesting point here is that as a consequence 
of the multiplicative noise in the amplitude equation, several state transitions may arise 
depending on the value of the noise strength. As it was pointed out in [31] , the presence of 
multiplicative noise, and in particular the fact that the coefficient o-^/j m A from Eq. (|2.6I) 
vanishes at A = 0, becomes crucial for the description of the asymptotic stability of the 
stationary PDF in terms of Lyapunov functions, giving rise then to different scenarios 
depending on the integrability of the PDF, i.e. whether the PDF can be normalised 
or not. As we shall show now, different states can be observed by simply changing the 
strength of the noise. The results are presented separately for each case of boundary 
conditions. 

In the case of DBC, the PDF of the dominant mode amplitude, A(t) = \u\(t)\, is 
given by Eq. ()2.9[) with A = 0. In this case there exist different states for the fluctuating 
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Figure 2. (Colour online) Numerical results for the noisy KS equation (|3.1|l integrated 
on a [— 7r,7r] domain with DBC (ID null space), (a) PDF of the first-mode amplitude 
calculated for a — 12, a — 34, and a = 80, with e = 0.05. Dashed lines correspond 
to a data fit to a function of the form given by Eq. (|2.9|) with A = 0. The inset shows 
the typical fluctuations of the amplitude at each of the three values of a. (b) Maximum 
location of the amplitude PDF at different values of a. The solid line corresponds to the 
numerical solution of the SL model. All values have been normalised to the corresponding 
value at o = 10. 



amplitude that can be characterised by looking both at the maximum location and 
integrability of the amplitude PDF [16j [31]. First, we observe that as long as ot\ > 
the maximum of P{A) occurs at a finite value, A mllx > 0, and then the state of A is 
characterised by finite fluctuations around a mean value (state I) . On the other hand, for 
— 1 < ot\ < 0, the maximum is located at zero, j4 max = 0, and the amplitude fluctuates 
intermittently between zero and a finite value (state II, see also below in §3.1.11 for a 
detailed statistical analysis of this state). These two states are separated by the critical 
value: 

(7i = ( 7m /2- 7 ir 1/2 - (3.4) 

Note that for 71 > 0, this transition can only be observed as long as "f rn > 271, while it 
is always observed for 71 < 0. By computing the maximum of P(A) at different values 
of a we can then characterise the critical behaviour as: A max = \a\ — a 2 ] 1 / 2 /(cri^/72) 
for a < <7i, and A max = otherwise, so that A max and a 2 are the order and control 
parameter, respectively, describing the critical transition. By using now the particular 
values obtained for the KS equation [from Eq. (|3.3|) we get 71 = —1/2688, 72 = 1/48, 
and 7 m = 1/576] we find <j\ — 28.3 in excellent agreement with the numerical results 
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Figure 3. (Colour online) Time averaged solution of the noisy KS equation solved with 
DBC for e = 0.1 and a — 10 (a) and 55 (b) for the case of: highly degenerate noise acting 
on the second mode (case I, solid line), degenerate noise acting on both the second and 
third mode (case II, dashed line), and highly degenerate noise acting on the third mode 
(case III, dot-dashed line). 



(presented in Fig. [2j where we obtain a\ ~ 29). Finally, we note that since 71 < 0, a 
second transition occurs when u\ < — 1. In this situation, it is not possible to define 
Lyapunov functions to show the existence of a globally stable stationary PDF for A 
which can be normalised [31] , As a consequence, the PDF converges to a Dirac delta 
function, P(A) = 5(A), describing then a completely stabilised state with A = (state 
III). The critical value an for this second transition corresponds to: 

trn = n/VW, (3-5) 

which for the KS equation gives an = 51.8, in very good agreement with the numerical 
results [see the bottom panel in the inlet of Fig. Ufa)]. Interestingly, such stabilisation 
effect is reflected on the zero solution obtained when the field u is averaged in time, 
(u(x,t)) t — [see solid line in Fig.[3fb)]. Since the dominant mode has been completely 
stabilised, the dynamics of the solution is given by the noise coming from the stable 
modes, which when averaged in time, gives zero solution. This is in contrast to the 
inhomogeneous solution observed for smaller values of a corresponding to state I [see 
solid line in Fig. Ufa)]. 

When we consider the case of PBC, the amplitude PDF is given by 

P(A) oc A a2 exp (-n'A 2 ), (3.6) 

where now ol-i = (1 + 27 1 <r 2 )/(27 m ) and /1' = 72/47 m . We first note that the transitions 
between different states can only occur if and only if 71 < 0. In such a case, the critical 
values for the first and second transitions are found to be: 

aj = VyWnl on = [2(|7i| - lm)}- 1/2 . (3.7) 

Noteworthy is that the second transition can only occur as long as 7 m < |7i|. Otherwise, 
the completely stabilised state III is never observed, and the PDF tends to P(A) ~ A a °° 
as a — y 00, with = — |7i|/7 m . By using the KS coefficients given in Eqs. Q3.3p . we 
obtain the first transition to occur at a\ — 36.3, in excellent agreement with the numerical 
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Figure 4. (Colour online) Numerical results for the noisy KS equation ()3.1[) integrated 
on a [— 7r,7r] domain with PBC (2D null space), (a) PDF of the first-mode amplitude 
calculated for a = 12, a = 38, and a = 80, with e = 0.025. Dashed lines correspond 
to a data fit to a function of the form given by Eq. (|3.6|) . The inset shows the typical 
fluctuations of the amplitude at each of the three values of a. (b) Exponent oli obtained 
from the data fit at different values of a. The solid line corresponds to the analytical 
solution of the SL model given by Eq. Q3.6[) . 



results presented in Fig. |4] The second transition, however, cannot be observed since the 
condition j m < does not hold, giving an asymptotic behaviour with aoo — —0.21, 
again in very good agreement with the numerical data [cf. Fig. [4]^6)] . 



3.1.1 Universal intermittent behaviour 

Let us focus here on the intermittent behaviour observed between the two transitions, 
g\ < a < a\\, for DBC. It is worth emphasising at this point that on-off intermittency is 
a crucial universal feature of many nonlinear systems close to criticality, and reflects a 
transition from order/coherence to a disordered state (hence understanding the statistical 
properties of intermittency is crucial for the characterisation of this transition) . In our 
case, it reflects the transition between an initially inhomogeneous state in space to a final 
zero state (see solid line in Fig. [3]). 

Figure [5] depicts the dynamics of the amplitude calculated by using a = 48, which 
for DBC is close to the second critical transition an. In this regime, fluctuations are 
clearly dominated by an on-off intermittent, or burst-like behaviour. As was pointed out 
in Refs. [15, 39 , such kind of on-off intermittency can be characterised by studying the 



Additive noise effects in active nonlinear spatially extended systems 



13 




t 

Figure 5. (Colour online) Time evolution of the dominant mode amplitude A(t) for DBC 
with a = 48. The waiting time between two consecutive large events is denoted as T. 
The dashed line represents the position of the chosen threshold Cth which defines the zero 
state as A < c t h- 




Figure 6. (Colour online) Statistical analysis of the waiting times between two consecu- 
tive bursts observed in the intermittent state II for DBC. (a) PDF of the waiting times 
calculated by using different noise strengths. The dashed line corresponds to a data fit 
to P(T) = NT- T exp{~T/T ) with r = 1.5. The value of the threshold used to define 
the waiting times T is Cth = 1. (b) Value of the fitted time scale To as a function of a. 
The solid line corresponds to the analytical expression given in Eq. (|3.9[) . All values have 
been normalised to the corresponding value at a = 12. 



statistical properties of the waiting time between two consecutive bursts, defined as large 
fluctuations above a given small threshold, A > cth- Given the SL amplitude equation 
(|2.6|) . we can obtain an analytic expression for the PDF of the waiting times, P(T), as 
follows. 

Let us assume an initial value which is below the threshold, i.e., A* = A(t = 0) < Cth- 
We then seek the probability P(T) that at time T, the variable A reaches the threshold 
Cth for the first time. In this close-to-zero state, the amplitude is found numerically to be 
A < 0.1, and therefore small enough to neglect the nonlinear term in the SL equation. 
More precisely, we first introduce the transformation y = log A, and assume that y <C 
to obtain the following linearised equation: 



y = K + cr^/^W(t), 



(3.8) 
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defined in the semi-infinite domain y G (—00, yo] with adsorbing boundary conditions 
at the origin yo = logcth, and where k = 1 + 71 cr 2 . We note that the choice of an 
adsorbing boundary condition at the origin is required to ensure that the variable A 
leaves the domain once it reaches the threshold. The problem now reduces to finding the 
probability that at time T the new variable y reaches the origin for the first time. This 
corresponds to the well-known "first-passage probability" (FPP) of the random walk in 
semi- infinite domains (see Ref. [H]), which simply reduces to solving the Fokker-Planck 
equation for the probability p(y, t): 

d t p(y, t) = -Kd y p(y, t) + ^ m <r 2 dyp(y, t). 

This can readily be done by using the method of images, obtaining: 

p(y,t;y l: 0) = 1 [p-fa-^+«*) 2 /4 7mCT 2 t _ ^/(^vf e -( w+l «+»rf) a /i Jm <r 2 t^ 
v / ir4:'y m cr 2 t 

where yi = log A4 is the initial value and we have assumed without loss of generality that 
cth = 1, and therefore yo = 0. The FPP can now be obtained as: 



P{T) = - [ d t p(y,t;y t ,0)\ t=T dy, 



from which we get: 



P(T) 



-(-y t +aT) 2 /4f„ 



\J 7r47 m cr 2 T 3 

In the long-time limit T — > 00, the above expression becomes: 

P(T) « T- 3/2 e~ T/To : (3.9) 

with a time scale To = [2a ^/j m / (1 +71CT 2 )] 2 . Therefore, at the critical point a = an the 
above PDF becomes a pure power-law dictated by a universal exponent 3/2 that does not 
depend on the particular model we are using, i.e., it does not depend on the coefficients 
of the SL equation. The numerical results obtained with the KS equation are presented 
in Fig. [6] where we can see an excellent agreement with the above expression. It should be 
noted that far from the critical point (cr < en), the power-law is exponentially corrected 
with the time scale T [cf. Fig. [6^6)]. 



3.2 Case II. Amplitude equation with both additive and multiplicative noise 

Let us now consider the noise term acting on both the second and third mode, i.e. — 
Sk.2 + <5/c,3 with DBC. In this setting, the amplitude equation is given by Eq. (|2.6|) with 
both 7 a and 7 m being different from zero. Although the noise is still highly degenerate, 
it induces the same effect as when we consider the noise acting on the whole subspace 
of stable modes, i.e., qu = 1 for k > 2 (and therefore having space-time white noise). 
For simplicity, we shall therefore restrict our analysis to the case with only the third and 
second modes being stochastically forced. 

First, we note that in this case of noise we can still identify two different dynamic 
behaviours (states) for the dominant amplitude. At low values of a we have that, as 
it occurred in case I, the location of the maximum of P(A) is found to be at a finite 
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Figure 7. (Colour online) Numerical results for the noisy KS equation (|3.1[) integrated 
on a [—it, 7r] domain with DBC and with a noise acting on the second and third stable 
modes (case II). (a) PDF of the first-mode amplitude calculated for a = 12, a = 36, and 
g = 80, with e = 0.1. The inset shows the typical fluctuations of the amplitude at each 
of the three values of a. (b) Exponent ot\ obtained from the PDF at different values of a. 
The solid line corresponds to the analytical solution of the SL model given by Eq. (|2.10[) . 



value, A max > 0, and the dominant mode fluctuates around a finite mean value. As the 
noise strength is increased, the maximum location approaches zero and the first mode 
component u\ may reach zero. Computing the maximum location by considering the 
PDF given by Eq. (|2.6p we obtain that the critical value, <j\, separating these two states 
corresponds to the same value as in case I of noise given by Eq. (|3.4j) . In this situation, 
however, the presence of the additive noise keeps the first mode from remaining at the 
zero position, and hence ruling out the intermittent dynamics observed in the previous 
section. Further, as the noise strength is increased, the PDF can now be always normalised 
so that the completely stabilised state III is never observed. 

Figure [7| depicts typical numerical results for the KS equation. The coefficients of the 
SL equation in this case are 7i = -1/2688+ 1/52416, 72 = 1/48, 7 a = 1/580608, and 
7, n = 1/576, so that we have o~\ ~ 28.6. We compute the amplitude of the dominant 
mode A{t) = for different strengths of the noise. The PDFs of A for a = 12, 

36 and 80 are presented in Fig. [71(a). We can see that the presence of the extra term 
controlled by 7 a in the amplitude equation keeps the dominant mode away from zero 
and hence of being completely stabilised. As before, we find excellent agreement between 
the numerical results and the analytical derivation given by Eq. (|2.6p . Indeed, Fig. [7j6) 
shows a comparison of the numerical value of the exponent ol\ obtained from a data fit 
of the PDF with the analytical value given by Eq. (|2.10[) for different strength of the 
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noise. Moreover, by computing the time average of the spatio-temporal solution we can 
see that at low values of the noise strength (a = 10), the averaged solution still retains 
the sinusoidal shape [see dashed line in Fig. QJa)] , corresponding to the first state where 
A max > 0. At higher noise strength values (ex = 55), however, the averaged solution losses 
the sinusoidal shape [see dashed line in Fig.[]T6)], corresponding to the second state with 
A max — 0. In contrast to the stabilised state observed in case I of noise, the averaged 
solution remains now noisy since the dominant mode is never completely stabilised. 



Finally, we consider the case of DBC with a noise acting only on the third mode (qk — 
S^k)- In this setting, both noise coefficients are zero (pf m = j a = 0) with a positive 
coefficient 71 > for any value of the disorder strength. As a result, the dynamics of the 
dominant mode is not affected by the noise, and only its steady state A st , given by the 
stationary solution of the deterministic SL Eq. (j2.6[) as: 



will be affected for sufficiently large a. When we apply this noise setting to the KS 
equation we obtain 71 = 1/52416, so that only for values of noise strength as large as 
a ~ 230 such effect will start to be relevant. As we see in Fig. [3] (see dot-dashed line), the 
time averaged solution (u(x,t)) is not affected when the noise is increased up to a = 55. 
It is also worth emphasising that we would obtain the same result had we considered a 
noise acting on all the stable odd modes. It means that this kind of noise acting on the 
odd modes only is filtered out by the quadratic nonlinearity interaction, such that it has 
no effect on the dominant mode dynamics. 



Let us consider now the situation where there are more than one unstable modes. Clearly, 
the theory presented in ^2 . II is not expected to be valid in this case. We shall therefore 
numerically integrate the noisy KS equation by taking e = 0.87 such that the number 
of unstable modes is two. We are interested to see whether stabilisation phenomena as 
presented in EI3.1l for noise case I can still be observed by appropriately tuning the noise. 
To this end we choose the noise so that it acts on the first stable mode, corresponding 
to taking qj~ — <5fc,3, and we study the dynamics of the amplitude of both the first 
(Ax(t) — \ui(t)\) and second (A 2 (t) = \u 2 (t)\) modes. 

The corresponding numerical results are presented in Fig. |8l Noteworthy is that both 
modes can be completely stabilised as the noise strength is increased up to values of 
a ~ 270, going through the same critical transitions as in 8 33.11 As before, we perform 
a statistical analysis of the waiting times between two consecutive bursts which are 
observed in the intermittent state corresponding to a — 260. The results for both modes 
are presented in Fig. HI where we observe that the PDF of the waiting times for both 
modes is dominated by a heavy-tail function with an exponent 3/2. Therefore, our results 
show that even in situations relatively far from the instability onset, we can still observe 



3.3 Case III. Amplitude equation without noise 




(3.10) 



4 Numerical results in a regime beyond the instability onset 
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Figure 8. (Colour online) Numerical results for the noisy KS equation (|3.ip integrated 
on a [— 7r,7r] domain in a regime far from the stability onset (e = 0.87) with a highly 
degenerate noise acting on the first stable mode. Panels (a), (b), and (c) show the time 
signal of the amplitude of the first dominant mode A±(t) for a — 40, 230 and 270, 
respectively. Panels (d), (e), and (f), show the corresponding time series of the second 
unstable mode A2(t) for the same noise strengths. Bottom panels (g), (h), and (i) show 
the PDF of each mode for the three different values of a. 




login CT) 

Figure 9. (Colour online) PDF of the waiting times between two consecutive bursts 
obtained from the time series of the two dominant modes A\ and Ai for a = 260. The 
dashed line corresponds to a data fit to P(T) = NT~ T exp (— T/Tq), with r = 1.5. 



stabilisation and critical transition induced by pure additive noise. We note, however, 
that the noise strength needs to be much larger now (up to a — 270) to completely 
stabilise the dominant modes. 
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5 Conclusions 

We have undertaken an analytical and numerical study on the effect of additive noise on 
active nonlinear SES. For simplicity we have focused on SES with quadratic nonlineari- 
ties, in particular the KS equation, and randomly perturbed in the vicinity of a primary 
bifurcation point, where the dynamics is described by a single unstable dominant mode. 

By adding a stochastic forcing acting on the first stable mode, we have been able to 
provide a detailed and systematic investigation of the transitions between different states. 
In particular, we have observed that at low values of the noise strength, the amplitude 
of the dominant mode is dominated by finite fluctuations around the zero-noise solution 
(cr = 0). On the other hand, at high values of the noise strength the dominant mode can 
be completely stabilised, yielding the zero solution when the spatio-temporal evolution 
of the system is time averaged. These two states are continuously connected by another 
intermediate state where the solution intermittently fluctuates between the first non- 
zero state and the final zero state, and it is characterised by a burst-like dynamics. 
The transitions between the different states have been completely characterised through 
critical exponents, obtaining excellent agreement between theory and numerical results. 
In addition, we have been able to rigorously derive the critical exponent 3/2 describing 
the waiting times of the intermittent state II. It is important to emphasise that the 
derivation has been done for general SES whose dominant mode amplitude is given 
by Eq. (|2.6j) . In this sense, the exponent 3/2 reveals the existence of an underlying 
universal mechanism dictated by the random walk properties. Interestingly, the same 
exponent has been ubiquitously found in many physical systems that display avalanche 
or burst-like dynamics. Examples include neuronal avalanches in the spontaneous cortex 
activity [46], on-off intermittency in electroconvection of nematic liquid crystals [21], or 
interface dynamics in disordered media [28] [40] . 

When the noise acts on the whole subspace of stable modes (or equivalently on the first 
and second stable modes) the dynamics of the dominant mode is corrected by an additive 
noise which keeps the system away of being completely stabilised, with the loss of any 
critical intermittent behaviour. Again, the numerical results in this situation are found 
to be in very good agreement with the analytical predictions. We have also considered a 
case where the noise acts on the second stable mode only (or equivalently on the stable 
odd modes). In this situation the quadratic nonlinearity can filter out the noise such that 
the dynamics of the dominant mode is not affected by the stochastic forcing applied to 
the system. 

It is important to remark that, although the theory presented in this work has been 
applied to the KS equation, our analysis is rather general and could be easily extended 
to other models, such as the stochastic Burgers equation [5j [43] , used for example as a 
prototype for ID turbulence albeit without pressure gradient, the Kardar-Parisi-Zhang 
(KPZ) equation, largely studied in the context of surface growth p], or the stochastic 
Swift-Hohenberg equation, often used as a model for Rayleigh-Benard convection, which 
in turn is commonly used as a prototype to study instabilities out of equilibrium in 
SES [7]. It should also be noted that the approach followed in this study is very robust: 
The stochastic SL equation (|2.6|) is the universal amplitude equation for SPDEs with 
quadratic nonlinearities, as proved in Ref. [4]. Furthermore, the characterisation of the 
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different regimes depending on the strength of the noise, and in particular the calculation 
of the critical exponent that describes the intermittent behaviour are systematic and 
rigorous. It is expected that the techniques and methodologies developed in this study 
can be extended to more general classes of problems, where noise induced phenomena 
occur. 

Finally, by performing numerical integrations of the noisy KS equation in a regime 
relatively far from the instability onset, where there are two unstable dominant modes, 
we have observed similar stabilisation and state transitions induced by an additive noise 
which is acting on the first stable mode only. It is important to remark that, although 
there is no theory in this case, our numerical study provides evidence that such stabil- 
isation and noise induced state transitions are not restricted to a regime close to the 
instability onset. 



Appendix A Derivation of the amplitude equation 

We use techniques from homogenisation theory [38j and singular perturbation theory for 
Markov processes [37] (see also [32]) to derive the amplitude equation that describes the 
dynamics of Eq. (|2.4|) near the bifurcation point. 

We consider a finite-dimensional truncation of Eqs. ()2.5[) up to N modes with a finite- 
dimensional null space, i.e. dim (J\f) = Nq, so that the number of fast modes is M = 
N — Nq. The kernel JV of C is then spanned by the first No cigenfunctions, and we can 
write: 

N N 

vi = /l Qfcefc and v± = ^ Vk^k, 

k=l k=N +l 

which when introduced into Eqs. (|2.5p yields the following system of equations: 

a m = -J% l {a,y)+f?{a), m = l,...N , (Ala) 

yk = -^kVk + -9o(a,y)+9i{y) + -(rqk$k, k = l,...M (Alb) 

where a = (ai, . . . , a./v ) T and y = (y\, . . . , 2/m) T , and for simplicity we have reordered 
the fast modes so that the label goes as k = 1, . . . , M. The different vector field terms in 
the above system of equations are given as: 

M M 

/o"(a, y) = 2a m B m e m ye + ^ B n £ m y n yi, (A 2 a) 

i=\ n,l=l 
N N M M 

go(a,y)= ^ B nik a n a e + 2 a n ^ B nik ye + ^ B nlk y n y e , (A 2 b) 

n,£=l n=l t=l n,e=l 

/i TO (o) = j m a m , g k i(y) = jkVk, (A 2 c) 

where the coefficients j k and are defined as Je k = jk^k, and — Ce k — \k^k, respec- 
tively, and 

B n ek = (B(e n ,e e ),e k ) , (A3) 
with the inner product (f,g) = J*" f(x)g(x) dx. We shall assume here the condition 
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V c B(ck, Ck) = for the bilinear map, so that we have 

B kkm =0, m=l,...,N . (A 4) 

This assumption, which is satisfied for the Burgers nonlinearity ud x u, ensures that the 
centering condition from homogenisation theory, see Eq. (|A 81) below, is satisfied. The 
system of equations (|A f [) for a m and y k is of the form of a fast-slow system of stochastic 
differential equations for which the associated backward-Kolmogorov equation (for w t = 
E(f(a(t),y(t)\a(0) = a,y(0)=y)) reads: 



d t w< = (e- 2 £ Q + e~ 1 C 1 + C 2 )w £ , (A 5) 



where 



A I 



fc=i ^ 

N M 

A = E /o m («> v)dm + E so(«. y)^, 

m — 1 fe= 1 

c 2 = E/r(«)^ + E-9o(y)^- 

m— 1 A;— 1 

and d m and 9fc represent derivatives respect to a m and j/a; for m = f , . . . , iVo and k = 
1,...,M, respectively. Here, the operator £q corresponds to the generator of a finite- 
dimensional Ornstein-Uhlenbeck (OU) process, so that the invariant measure of the fast 
process is Gaussian: 

f M X 

p(dy) = -? ex p(-E -r^yljdy, ( A7 ) 

where Z is the normalisation constant. From (IA 2 a|) and Assumption (|A4[) we deduce 
that the vector field /™(a,y) is centered with respect to the invariant measure of the 
fast process, 



fZ l (a,y)p(dy) = 0. (A 8) 

We will now show that the amplitude equation is of the form 

dz = v M (z)dt + g M (z)dW, (A 9) 

where W denotes a standard A^ — dimensional Wiener process. The subscript M used 
here is to emphasise the fact that the homogenised coefficients depend on the number of 
fast processes that we take into account. The calculation of the coefficients vm(z) and 
g M (z) which appear in Eq. (|A 9p will be obtained by using homogenisation theory. 

We solve the backward-Kolmogorov equation (|A 5P by looking for a solution in the 
form of a power series expansion in e: 

w e =w + ew 1 +e 2 w 2 +O{e 3 ). (A 10) 

Substituting this expansion into (|A 5j) and equating coefficients of powers of e to zero, 
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we find: 

£ w =Q, (Alia) 

£o^i +£iw = 0, (A 11 6) 

£oW2 + C\w\ + C2WQ — d t WQ. (A 11c) 

Since the null space of the generator of the OU process consists of constant in y, from 
Eq. (| A 1 1 a[) we deduce that wq — wo(a,t). The solvability condition for the second 
equation reads 

/ (£ lWo )p(dy) = 0. (A 12) 

JR M 



This solvability condition is satisfied on account of (I A 8[) . and the fact that the first term 
in the expansion is independent of y. We now solve this equation by using the following 
ansatz: 

, M M s 

wi= {a m '^2xe m ye+ ^ (ne m ynye)d m w , (A 13) 

^ e=i n ,e=i ' 

where the constants xe-m and (, n im are to be determined from Eq. (|A 11 6[) together with 
(|A 12|) . In particular, we obtain: 

2-Sm/cm u ^nkm ( K \ A\ 

Xkm = \ ! Ufcm = T —5—) (, A14 J 

for m = 1, . . . , No and n, k = 1, . . . , M. Finally, for Eq. (I A 11 cl) we have the solvability 
condition: 

/ ( d t w - £ 2 wo - A^i ) p{dy) = 0, (A 15) 

JR M \ J 

which allows us to obtain the homogenised SDE (|A 9lFI . For the components of the drift 
term we get: 

/ , M M M \ , 

»M («) = (( fi + a m J2 9oXkm + 2 doVnCnkm + /<T X fe ™^ ) ) ' ( A 16 ) 

^ ^ k=l n,k=l k=l ' ' 

for m = 1, . . . , iVo, and where (...) denotes average with respect to the invariant measure 
p(dy). Similarly, the components of the quadratic form associated with the diffusion 
matrix g M are given by: 



2 (Sm)« = (fa [ a j Xn i Vn + CntjVnVe] ) , (A 17) 



where i,j = 1, . . . , Nq. The integrals in (|A 16|) and ([A 17|) are Gaussian integrals that can 
be calculated analytically. In particular, substitution of /g and g l given in Eq. (|A 2 a[) 
and (IA 2 respectively, into (|A 16|> and (|A 17|) . and calculation of the resulting Gaussian 
integrals lead to the amplitude equation (IA 9[) which can be written as the following 



4 To be more precise, it allows us to obtain the homogenised backward Kolmogorov equation 
from which we can read off the limiting SDE. 
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system of coupled SL equations: 

N 

da m = VjJHai, . . . , a No )dt + ^(<7Ar) mn <iW n , (A 18) 

71=1 

for m = 1, . . . , -/Vq, and where the different terms are given as: 



M od r> \ ^0 M OR R 

V <-'- D kmni- D nnk 2 2 \ , \ ^ \ T i -'- D nkm- D £nk o 2 



JV M OR R 

E\ ^ ^-Okmm- D n£k , . s 

a„o; > , (A 19 a) 

Ai- 

n,t=l k=l K 

1 M 1R R M R R 

/-2 "\ ^- D lii JD jnn 2 2 i JD kni JD knj 422 /A in l\ 

^Wm^ = a^i 2^ CT 9£ + \, 7T I X^2 g 9n9 fe - (A 19 b) 



In the case of a one-dimensional null space (No = 1) we define A = a%, and the above 
expressions have a simpler form. In particular we have: 

v M (A)^( J1 +- fl a 2 )A~ l2 A 3 , (A20a) 

g M {A) = <y^J la a 2 + lm A\ (A 20 b) 

with the different coefficients given as: 

N N N 

E2B n iiB n ik 2 , B n nB nn £ 2 \ 2B n nBu n . 

{Xn + X e )A n q " + 22 A^Ai 72 = E A n ' (A21a) 

n,fc2 v 7 n,fc2 n=2 " 

N OR2 W AR2 

\^ LD nl\ „2 „2 „, V *^"11„2 /» 01 M 

„,fc2 ^" + Xl > Xi n=2 A « 

where the noise in the resulting amplitude equation is interpreted in the Ito sense, and we 
have reordered the fast modes back again so that the label goes from k = No + 1, . . . , N. 
It is also interesting to consider the case when the null space is two-dimensional. By 
assuming the Burgers nonlinearity (ud x u), we can explicitly work out the different terms 
of the amplitude equation (|A 181) for the two components ai (t) and a 2 (t) , obtaining: 

Vm = —TT a l ~ 7^- ai °2 + O'l + v M)a x , (A 22 a) 

4A2 4A2 

v 2 M = --pr-a! - T^afa 2 + (j 2 + V M )a 2 , (A 22 b) 

4A2 4A2 
/ ^ rial ^ q %a\ r a i &J± n n 

\<?m) = T "*% t . - I ' A 22 , , 



with the coefficients: 



Vm _a 2 ™ k\ k (gl +1 + rg +1 ) - \ k+1 (<ll + r|)(fc + 1) 
8 ~* 2 (Afe+i + A&)AfeAfe+i 



4 jvj 



<7fc<7fc+i + r Ffe+ 



16 ^-^ Afc(Afe + i + Afc)Afe + i ' 
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Figure B 1. Sketch of a liquid film falling down over an inclined plane with a vibrating 
disordered wall. The thickness of the film and the disordered wall position are denoted 
as h(x,t) and s(x.t), respectively. 

where we have assumed a general case where the two components of the noise for each fast 
mode (which arise as a consequence of the fact that the null space is multidimensional) 
have different wave number dependence, namely qk and r^. We therefore see that in 
general, the diffusion matrix is non diagonal and the resulting amplitude equations are 
not coupled only through the drift term but also through the noise term. 



Appendix B Derivation of the noisy Kuramoto Sivashinsky equation for a 

hydrodynamic system 

We provide a physical example where the noisy KS equation can be derived from a general 
formulation in the framework of thin-film hydrodynamics. 

Consider gravity-driven thin film flow down an uneven wall that is inclined at an angle 
9 to the horizontal, as shown in Fig. IB ll Let x and z be the stream- and the cross-stream 
coordinates, respectively, and the wall be described by z — s(x,t) relative to a datum 
plane and where the time-dependence is such that each point of the wall fluctuates only in 
the z-direction, i.e., at a given time t the velocity of the wall point [x, s(x, t)) is the vector 
(0, s t (x, t)). Let also the free surface be located at z = fix, t). We introduce dimensionless 
variables by utilising the Nusselt film thickness ho (corresponding to undisturbed flow 
down an inclined plane) as a lengthscale, the Nusselt surface velocity Uq = pgti^ sin 6/2(1, 
where p and (i are the density and the dynamic viscosity of the liquid, respectively, and 
g is gravity, as a velocity scale and Uo/ho as a time scale. The pressure scale is chosen 
as piUo/ho- We additionally introduce the Reynolds number Re = pUoho/p and the 
capillary number Ca = pUo/j where 7 denotes the surface tension. 
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The governing equations are given by the incompressible Navier-Stokes equations: 

Re(u t + uu x + vu z ) = —p x + u xx + u zz + 2, (B 1) 

Re(v t + uv x + vv z ) = -p z + v xx + v zz - 2 cot 9, (B 2) 

u x + v z = 0, (B 3) 

where subindices represent partial derivatives, u and v denote the x- and the z-component 
of the velocity in the liquid, respectively, and p denotes the deviation of the pressure from 
the atmospheric level. The boundary conditions at the wall, z = s(x, t), are given by: 

u = 0, v = s t (x,t), (B4) 

and at the free surface, z = f(x,t), the kinematic compatibility and the tangential and 
normal stress balance conditions are satisfied: 

ft + uf x -v = 0, (B 5) 

(1 - fl){u z + v x ) + 2f x {v z - u x ) - 0, (B 6) 

p = YTp x hfx[Uz + Vx) + Uxfx + Vz) " CaJiTflW^- (B 7) 

We consider here the long-wave approximation so that we introduce the so-called thin- 
film or long- wave parameter e, defined as the ratio of the typical film thickness to the 
length scale over which streamwise variations occur. We then introduce new variables: 
£ = ex, t = et, and w = v/s. We also assume that Ca = 0(e 2 ) and define Ca — Ca/e 2 . 
We expand the different variables as, u — Uq + eu\ + ■ ■ ■ , w — wq + ewi + ■ ■ ■ , and 
P — Po + £ Pi + ' ' ' 5 which lead to a series of solvable perturbation problems. The leading- 
order problem is: 

u zz = 0, poz = -2cot6», u o? +woz=0. (B8) 

At z = s(£, t), mo = 0, wq = s t (£,t). At the free surface, z = /(£, t), the tangential 
and normal stress balance conditions imply uq z = 0, po = —f^/Ca. The solution of the 
problem at leading order is: 

u Q = -(z-h-s) 2 + h 2 , (B9) 
w = u (h + s) s - [h\(z -s) + s r (£, t), (B 10) 

pa = -2{cot9)(z -h-s)- (h + s) K /Ca, (Bll) 

where we introduced for convenience the film thickness h = f—s. Note that the kinematic 
compatibility condition can be written as 

h T + q c = 0, (B 12) 
where q = J u dz. Using (|B 9p . we find q = |/i 3 + O(e), which implies 

h T +[lh 3 ] 6 + 0(e) =0. (B13) 
At next order, we obtain the following system of equations: 

u\zz = Re(u 0r + u Q u oi + w Q u Qz ) + p £, (B 14) 

Piz = w 0z , (B15) 

ui^+w lz ^0, (B16) 
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subject to tii = w\ — at z — s and u\ z = 0, p\ = 2wq z — 2f^u^ z at z = f, where uq, wq, 
and po are given by (|B 9[) - (|B lip . The time derivative uq t involves the time derivative 
h T , as is evident from (|B 9p . The latter is eliminated by using (|B 1 3|) . For brevity we 
do not show the solution of the problem at first order (which can be easily found with 
a symbolic manipulation software) as it turns out to be rather lengthy. Substituting 
u = uq + eui + 0(e 2 ) into q — J u dz, we find 



2 

q = -h 3 + eh 3 



8i?e l3L 2 cot0,, , 1 „ . 
— h 3 hs —(h + s) s + -^-(h + s) m 



Thus, we obtain the following evolution equation for the film thickness, h: 

= 0. (B17) 



!>-,- [ ^h 3 +eh 3 



8i?e L3L 2cot6». L . 1 , 
— h 3 hz —(h + .s) ? + -=-(h + s) c « 



We proceed next with a weakly nonlinear analysis. We assume that both the amplitude 
of the free surface and the bottom wall are small, of 0(e), and we write h = 1 + erj and 
s = ea. Substituting these expressions into the evolution equation above re-written in 
the moving frame, \ = £ — 2r, we obtain the following evolution equation: 

fjf + 4f?f? x + Dfj xx + -^-rjxxxx = s - ( B 18 ) 

where 0(e 2 ) terms have been neglected, f — er, fj(x, f) — r)(x + 2f/e,f/e), r) = 
a(x + 2f/e,f/e), D = 8Re/15 - 2cot6»/3 and S = (2 cot 8/3)a xx - [l/3Ca)a xxxx . We 
find that s(£, r) = ea(£ — 2t, er), which physically means that the typical deformation of 
the topography shape is of small amplitude and large wavelength (of the same order as 
the wavelength of the typical free-surface wave) , it propagates downstream with constant 
velocity, and is slowly changing in the frame moving with this velocity. 

To simplify the latter evolution equation, we introduce the transformation, fj = U/A, 
X = X/B, and f = T/C, where A = 4/ (3Ca\D\ 3 )^ 2 , B = (ZCa^) 1 ' 2 , and C = ZCaD 2 , 
which leads to the following canonical form: 

U T + UU X ± U XX + Uxxxx = S(X, T), (B 19) 

where S — (A/C)T,, and the sign +/— corresponds to the positive/negative value of D. 
We therefore obtain the noisy KS equation (|3.ip . 
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